Spin excitonic and diffusive modes in superfluid Fermi liquids 
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A role of a particle-particle p-wave spin interaction in Fermi liquids with s-wave pairing is studied. 
Depending on the sign of the interaction there arises either the new exciton collective mode below 
the pair-breaking threshold or the diffusive excitation mode above the threshold. The Landau 
parameters which control the interaction strength are evaluated for various systems: the dilute 
fermion gases, degenerate electron liquid, metals, atomic nuclei and neutron matter. The interaction 
removes also the square-root singularity in the phase space of pair breaking processes. It is shown 
how these effects influence the neutrino emissivity in the neutron Cooper-pair recombinations in 
neutron stars. 



Processes with recombinations of Cooper pairs pro- 
vide important information about interparticle interac- 
tions and the pairing mechanisms in different fermionic 
systems: ordinary superconductors [l|, liquid ^He and 
•^He-^He mixtures jil, cold atomic gases [1] , atomic nuclei 
neutron stars and other systems. In supercon- 

ductors, they are studied by absorption of infrared radi- 
ation or by the Raman scattering [7[ . In the cold fermion 
atom gas one can use the Stokes scattering to detect the 
onset of the pairing fS\ . Inverse pair breaking and forma- 
tion (PBF) reactions constitute an important mechanism 
of the neutron star cooling [1, 0| . In these processes the 
energy is released in the form of neutrino-antineutrino 
pairs radiating off the star. Superburst ignition depth is 
sensitive to the value of the PBF emissivity in the in- 
ner neutron star crust In 1I| the PBF processes 
are suggested to be responsible for the recently observed 
rapid cooling of the young neutron star in Cassiopeia A. 

It was shown 12| that a residual interaction of single 
particle excitations, which does not contribute to pair- 
ing, can bind them in a state orthogonal to the Cooper 
pair, generating collective excitation modes in supercon- 
ductors. For superfluid ^He the similar mechanism was 
studied in [13|. Interactions in the same spin channel, in 



which the pairing occurs, were studied so far. The influ- 
ence of the interaction in one spin channel on the pairing 
in another channel has not yet been considered. 

In this Letter we study the effects of the p-wave inter- 
action in the spin-one channel on excitations in a Fermi 
system with the spin-zero pairing. We calculate response 
induced by the external spin- and helicity-density sources 
and show that depending on the sign of an effective in- 
teraction there appears either a new exciton mode or a 
diffusive excitation mode. Then we evaluate the strength 
of this effective interaction for different Fermi systems 
and, as an example, calculate the neutrino emissivity in 
the PBF processes for the neutron star with the neutron 
pairing in ^Sq state taking into account the effects of the 
new collective modes and correlations. 

We use the Fermi-liquid theory approach extended 
to systems with pairing by Larkin and Migdal and by 



Leggett [IJ]. For the processes induced by weak nucleon 
interactions this approach was adopted in Ref. [l^. In- 
teractions in particle-particle (^) and particle-hole (cj) 
channels are essentially different. The interaction am- 
plitude of two fermions with momenta p = pf n and 
p' ~ ppn' before and after the interaction in the 
channel is parameterized as = rg(ri, n')(itT2)(ic2) + 
Tf{n,n'){ia'2a'){aia2) and in the w-channel as F"^ = 
r^(n,n')l'l + ^i{n,n')a' a . Here pp stands for the 
Fermi momentum, n and ft' are the unit vectors. The 
unit matrices 1 and 1' and the Pauli matrices a and 
(?' act in the nucleon spin space. Superscript "w" in- 
dicates that the amplitude in this channel is taken for 
|(fw f| ^ w ^ ep, wf is the Fermi velocity, ep is the Fermi 
energy and q ~ {^,0} is the transferred 4- momentum. 
The coefficients of harmonic expansion of the scalar Fq'"^ 
and spin F^'" amplitudes, the Landau parameters, should 
be either evaluated microscopically or extracted from 
analysis of the experimental data 

The singlet pairing in the Fermi liquid occurs owing 
to the attractive interaction, a^pFg — Jq < 0. At 
zero temperature the paring gap A follows from equa- 
tion — l//o = Ao/{o-'^ p) ~ hi(2ep/A), where a is the 
pole residue and p is the density of states at the Fermi 
surface. This expression is naturally generalized for fi- 
nite temperature T, cf. Eq. (5) in the second paper of 
Ref. 15 1 . Since = for the scattering of identical 
fermions, the spin interaction in the ^-channel simplifies 
as a'^pT\{n, n') = g\(nn'). It is usually assumed that the 
higher Legendre harmonics are much smaller jl^. Since 
we focus on the spin channel, the interaction Fq decou- 
ples and can be dropped. In the cj-channel, a? pVf — 
9o ~^ 9ii^'^)- Contributions from the zeroth harmonics, 
.gg , are accompanied by the factor Vp, see Ref. p^], and 
for non-relativistic Fermi liquids under consideration (for 
Wp ^ 1) can be put zero. Thus we remain with only tree 
relevant Landau parameters /g < 0, and g". Let us 
first put gi zero and demonstrate the influence of the 
interaction in the ^-spin-one channel, g^, on the pairing 
effects in the ^-spin-zero channel. Then we recover de- 
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pendence on g". 

Consider now an external perturbation, which cou- 
ples to the spin density operator s{x) = ip^ {x)aipix) 
and the helicity density operator h{x) — ^\x){a'p + 
p<7)ip{x)/{2 m), where ijj is the spinor of a non-relativistic 
fermion, pis the momentum operator, and m is the mass 
of the free fermion. From these quantities one can build 
the axial (A) fermion current J"^ — {h,s). The Fourier 
transform of its bare components after the Fermi liq- 
uid renormalization becomes J'^''^{n,q) = i<^T^ i,<^t'^ q)- 
Here g = e^i/a and f^^ = CAVpn/a are w-vertices, 
CA is the effective charge of the quasiparticle. For 
F" = 0, that we now exploit, the in-medium vertices are 

TA,o{n,Q) = Ta,o and ta,i(^) = TajK?)- For F" 7^ 
these vertices are modified [15| . Additionally, in a system 
with pairing there arise new vertices responsible for the 
PBF processes: (jTA,o{n,q) and or A.i{n,q)- The y fo llow 
from the solution of the Larkin-Migdal equations |l^,[l5j, 

TAfl{ri,q) = -^({{nn'){N{n\q)+Ao)TAAn',q))r,. 

+ ((nn')0(n',g;-l)rXo)^,), (la) 
= ^^(^{{nn'){N{n',q)+Ao)¥AA{n',q))-, 
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+ ((nn')0(n',g;-|-l)fX,(n')),-j, 



(lb) 



The brackets indicate the angular averaging (. . . )fj = 
/ (...)• The loop functions 0{n, q; ±1) = ^a^p (z+± 
z_)gT{n, UJ, q), and N{n, q) — a^p z_ .9T(f?, q) with 
z± = (lo ± vq)/(2 A), and the master function 
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where E± ^ e+ ± e^, = /(e_) - /(e+), f{x) - 
l/(exp(a;/r) -h 1) and e± = [(^p ± vq f + A^jVs. The 
solution of Eq. ^ is 
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TAA = f n,?Si,o - ^7i(5T(n')i[l - (", n')'])fl'P± , 



where P± = ft — fiq (nftq), fiq ~ q/\q\ and the correlation 
factors 

bir' = iCo + (^^^^ip^5T(n)i[l-(nn,)2])_, 

[7«]-^ = iCo + (^^^^i^5T(n)(nn,)2)_ (3) 

are controlled by one effective interaction parameter 

Co = 3/g« - l/fl (4) 

The singlet pairing occurs for /| < and 3/q < g\. 
Then, if 5 J < 0, we have Co < 0, otherwise the p-wave 
pairing is preferable. For g\ > Q we have Co > 0. 



The response of the Fermi system to the excitation 
{A) is determined by the symmetrical current-current 
correlator n'^''(g) = i(Tr{ J'^''^(n, q) r{n, q)])n with the 
in-medium current J'^{n,q) = {^XA,i{n,q),(TXA,o{''^i<l)) 
expressed via the reduced current correlators derived 
in H: 

XAfii"^, q)^L{n, q; -1) TAfi{n, q) + M{n, q) TA,o{n, q), 
XA,i{n, q)^L{n, q; +1) taA^, q) + M(n, q) TA,i{n, q), 

where M{n,q) = -a^ p z+ gT{n,uj,q ), and ^'^'^J'f^^ = 
(l^ ''^)9T{n, {vcf),q)~{^- ^) grin, u, q) . The tem- 
poral and spatial components of the tensor are — 
{TA,iXA.i{n,q)}n, and = S'^TXnXAfi{n,q))n with 
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[grin, {vq),q) - gT{n,uj,q)'_ 
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The mixed components are 11'° — H"' = fiq J2j H"'"' • 

From Eq. ^ we see that the external perturbation can 
induce a singular response in the PBF amplitudes at the 
values ll) and q corresponding to the poles of the functions 
7^ and 7^ . These poles determine the new transverse and 
longitudinal collective modes (spin excitons). For 9 = 0, 
the longitudinal and transverse modes coincide and their 
frequency cj follows from the condition 
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0, y = c^/(2A), 



(6) 



where griy) = 5t(0, 2Ay - i 0, 0) . 

Although the full inclusion of the f;"-dependence is 
rather tedious, the modification of Eq. ^ is simply given 
by the replacement ^griv) ^griv)! (I + 5 5" ^9T{y)) ■ 
For I Co I ^ 1 it induces the shift 



Co^C = Co/(l + Co.9r/3). 



(7) 



This relation interpolates between the limits |Co| <C 
'i/\gX\ when C « Co, and |Co| > ?,/\gt\ when C ~ 
{?>/gi){l ~ 3/(g"Co)). So, parameter C controls effects 
of residual interactions on the PBF processes. 

In the long wave-length limit (for w > |g|) from Eq. © 

we get Sn*J (g) = ^5^3n00(w). The response function, 
having for y 
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FIG. 1: Left: energies of the collective modes (for g = 0) 
as functions of the parameter C for various values of A/T. 
Right: the response function R{y) for T = and y > 1 given 
by Eq. ^ for various values of the parameter C (solid lines). 
Dashed line shows the function R{y, C — >■ oo) . 



determines the probability of PBF processes. The cross 
section of the excitation scattering in matter is deter- 
mined by this response function R. 

Solutions of Eq. ^ are shown in Fig. [1] (left) as a 
function of C. For C < 0, solutions with y < I corre- 
spond to the undamped spin exciton branch at lo < 2A, 
since 55t(j/ < 1) = 0. For C > 0, solutions exist only if 
C > -K5t(1 + 0) and uj > 2A. Since here ^griv) + 0, 
they constitute the diffusive spin mode. The frequency 
of the exciton mode increases with the increase of T and 
decreases for the diffusive mode. The response function 
R(y) at T = and y > 1 is plotted in Fig. [T] (right panel) 
by solid curves for various values of the parameter C. For 
y > \ the function R{y > 1,C) only weakly depends on 
the sign of the value C therefore we present it only for 
C > 0. The function R{y,C oo) — "^gxiy) is shown 
by the dashed curve. It has a square-root divergence at 
y — ^ I + 0, which is smeared out for finite values of C. 
Thus the finite value of C leads to a reduction of the 
spin response of the Fermi liquid close to the threshold 
for a; > 2 A. A similar effect was discussed in Ref. for 
the Raman scattering on metallic superconductors. 

To estimate the value of our key parameter C we need 
to know Landau parameters gf, /g and g". For a di- 
lute Fermi gas we can use quasi-particle amplitudes de- 
rived in up to the second order in the parameter ^ — 
2aeffPF/Tr, where Oeff is the effective scattering length. 
We derive 4 = C+C^ (21n2-hl)/3, = 3C2(l-21n2)/5, 
= -2C2(ln2 + 2)/5 and obtain C « -5.7/(aeffPF)^. 
For the neutron gas the value of the nn-scattering length 
is very large, a ~ 20 fm, whereas the effective scattering 
length is much shorter |17| . being determined, e.g., by 
the Viow-'fc potential, UcS — 2 fm. 

For more complex systems the parameters in the ^- 
channel can be estimated with the help of the w-Landau 
parameters in the s~p approximation of 18] 
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FIG. 2: Parameters /g and gl (left panel) and C (right panel) 
for the neutron matter reconstructed with the tj-Landau pa- 
rameters calculated in Ref. [13, using Eq. ((9| as functions 
of the Fermi momentum. 



here Af 



/r/(i 



/r 



— Ji / ~^ jT+iJ — .y; / v-^ I 2iq 

The Fermi-liquid approach was applied to the degener- 
ate electron liquid in Ref. [l^ ■ Using Table I and Table II 
of [19] we find C — —2.54, e.g., for a small value of the 
parameter qbPf — 0.032, where ob is the Bohr radius. 

For alkali metals at zero pressure the first three w- 
harmonics are calculated in ReL 20]. Applying ^ we 
then find for sodium /|(Na) = -0.11, (Na) = -0.38 
and g"(Na) = —0.075. Here the p-wave paring is re- 
alized, since Co > 0, but the value \Co\ is very small. 
Bearing in mind large uncertainties in estimates of the 
w-Landau parameters one cannot exclude that C < at 
C| <g; 1. In the latter case we would deal with very pro- 
nounced effects of the spin exciton mode. This case can 
also be realized, if one allows a variation of the pressure. 
Thus presence or absence of the new exciton mode could 
tell about the kind of pairing in the given system. For 
potassium /|(K) = —0.56, (?i(K) = —0.89 and, using 
g^{K) = -0.12, we obtain C = -1.48. 

For the nucleon matter several harmonics of the lo- 
Landau parameters were evaluated in many works, e.g., 
see Refs. [l7, 21[. The parameter /q related to the ISq 
pairing was also calculated, see p. Contrary, the gf 
parameter is poorly known. Using results 17|, |21| we 
reconstruct g^ and /q with the help of Eqs. dH) and eval- 
uate then parameters Cq and C. For the neutron matter 
the results are shown in Fig. [2] in dependence of the Fermi 
momentum. We see that estimations of C are very un- 
certain due to discrepancy in different estimates of the lu- 
Landau parameters. Presented results show that might 
be |C| < 10-20 at some densities in the range of the ISq 
paring and even C might cross zero. Existence of regions 
where C < implies a possibility to observe effects of the 
exciton modes. 

Using the values of the w-Landau parameters and their 
density dependence extracted from the atomic nuclear 
experiments 0, [23]j we obtain C —10 for pF ~ Ifm"^. 
Thus the exciton mode could manifest itself in the nuclear 
surface phenomena. 

Now we apply Eq. ^ to calculate the neutrino emis- 
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FIG. 3: The ratio evvjsvl as a function of T/Tc for various 
values of C (see curve labels). 



may exist in superconducting potassium, in rare fermion 
gases, and in the neutron matter. In atomic nuclei the 
new spin exciton mode may manifest in the surface layer. 
Modification of the neutrino emissivity due to presence 
of spin excitonic and diffusive modes may have an impact 
on the neutron star cooling. 
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sivity in the neutron star matter in the region of ^Sq 
pairing. It is mainly determined by the neutron PBF 
process induced by the axial- vector current oc [l5| : 
the vector current contribution is 0(y\) and can be ne- 
glected [l^, [i^l • For one type of neutrino the emissivity 
then is given by 15 1 



'd\q\q^ iq^q,~g^,)^m-{q) 



48 n* 



exp(a;/r) - 1 



where G and qa are the weak-interaction and axial- vector 
coupling constants. The integration over \q \ yields 



35 7r3 



(10) 



where according to Eq. ^ there can be two contributions 
to e^p: one, for arbitrary C, from the pair-breaking con- 
tinuum with the diffusive modes at > 2 A and the other 
one, for negative C, from the spin-exciton mode with the 
frequency io{q) at < uj{q = 0) < 2 A. The later contri- 
bution is associated with the processes of breaking and 
formation of spin excitons. In the limit |C| — > oo the 
collective mode contribution vanishes as oc 1/|C| and we 
recover the result [l^ ; ; which follows from (|10p after 
the replacement R{y, C) — >■ R{y, G — >■ oo) = Qgriy)- 

Effect of the finite value of C on the neutrino emissiv- 
ity in the neutron PBF process is illustrated in Fig. [31 
where we plot the ratio Ei/p/eJ^p^ taking into account the 
standard temperature dependence of the ^Sq pairing gap 
A(T) ~ 3.1 Tc (1 - T/Tc)^/'^ with as the critical tem- 
perature. For \G\ - 5-10, cf. Fig. H (right), the effect 
becomes pronounced for T/Tc < 0.5, yielding a suppres- 
sion for C > and an enhancement for C < 0. Thus 
in different density regions there may arise either an en- 
hancement or a suppression of the PBF emissivity. Effect 
becomes even more pronounced for smaller values of |C|. 

In conclusion, we found that the spin p-wave inter- 
action in the particle-particle channel can produce new 
spin excitonic and diffusive modes in the Fermi system 
with the singlet paring. This interaction leads also to 
smearing out of the threshold singularity in the Cooper- 
pair breaking reactions. We calculated the relevant cou- 
pling parameters for several Fermi systems. Spin excitons 
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